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New measurements of Df and D* + meson production rates from B decays and from qq continuum 
events near the T(4S) resonance are presented. Using 20.8 fb~ of data on the T(4S) resonance 
and 2.6 ftT 1 off resonance, we find the inclusive branching fractions B(B — > DfX) — (10.93 ± 
0.19 ± 0.58 ± 2.73)% and B(B -> D*+ X) = (7.9 ± 0.8 ± 0.7 ± 2.0)%, where the first error is 
statistical, the second is systematic, and the third is due to the Df — > <f>ix + branching fraction 
uncertainty. The branching fractions EB(B -> d£" )+ .D ( * ) ) = (5.07 ± 0.14 ± 0.30 ± 1.27)% and 
EB(B -> D* s + D { * ) ) = (4.1 ± 0.2 ± 0.4 ± 1.0)% are determined from the D ( s * )+ momentum spectra. 
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PACS numbers: 13.25.Hw, 13.25.-k, 14.40.Nd 
I. INTRODUCTION 

The decay of B mesons into final states involving a 
(*)+ 

D s ' provides an opportunity to study the production 
mechanisms for cs quark pairs 1 . Although several di- 
agrams can lead to D s ' production in B decays, the 
dominant source is expected to be external W + — > cs 
emission (Fig. |l|). A precise knowledge of this produc- 
tion rate remains interesting in light of continuing the- 
oretical difficulties Q in accounting for the measure- 
ments of both the semileptonic branching fraction and 
the inclusive charm production rate in B decays. In- 
deed, it has been noted that an enhanced B decay rate 
to charm would help explain the small observed semilep- 
tonic rate 0. 

It is possible to produce D\ mesons in qq events 
from continuum e + e _ annihilation. The process of 
fragmentation (i.e., formation of hadrons) is non- 
perturbative and can only be modeled phenomenologi- 
cally. The ratio of vector to pseudoscalar production 
rates is of particular interest for testing such models. The 
Df system is well suited to measure this quantity because 
the cs states with L = 1 have not been observed to decay 
to either Df or D* + mesons. 

In this report, measurements of B — > DfX and B — > 
D* + X production rates and momentum spectra are pre- 
sented. We also determine the production cross section 
for Df and D* + mesons in continuum events. 
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FIG. 1: The main spectator diagram leading to the produc- 
tion of L>i*' + mesons in B decays. 



*Also with Universita di Perugia, Perugia, Italy 
^Also with Universita della Basilicata, Potenza, Italy 
[1] Reference in this paper to a specific decay channel or state also 

implies the charge-conjugate decay or state. The notation D$ 
means either or D$ + . B —* Dg ; s a general repre- 

sentation for any of the modes with cs and cq states including 
their excited states. The notation B — » Di*^ + X also implies 
B~ — » D { * )+ X. 



II. THE BABAR DETECTOR AND DATA SET 

The data used for this analysis were collected with the 
BABAR detector Q at the PEP-II asymmetric-energy col- 
lider || at the Stanford Linear Accelerator Center. An 
integrated luminosity of 20.8 fb _1 was recorded in 1999 
and 2000 at the T(45) resonance ( "on-resonance" ) cor- 
responding to about 22.7xl0 6 produced BB pairs, and 
2.6 fb" 1 at an energy about 40MeV below the T(45) 
mass ("off-resonance"). A detailed description of the 
BABAR detector can be found in Ref. Q . Only the com- 
ponents of the detector most crucial to this analysis are 
summarized below. 

A five-layer double-sided silicon vertex tracker (SVT) 
and a 40-layer central drift chamber (DCH) filled with 
helium-based gas are used to measure the momenta of 
charged particles. The tracking system covers 92% of the 
solid angle in the center-of-mass frame and lies within 
a 1.5-T solenoidal magnetic field. For charged-particle 
identification, ionization-energy loss (dE/dx) in the DCH 
and SVT, and Cherenkov radiation detected in a ring- 
imaging device (DIRC) are used. Photons are identified 
and measured by a CsI(Tl) electromagnetic calorimeter. 



III. THE D+ AND D*+ SELECTION 

Only the decay mode Df — > </>7r + with <p — > K + K~ 
is used since it has the best signal-to-background ratio. 
Charged tracks are required to originate within ±10 cm 
of the interaction point along the beam direction and 
±1.5 cm in the transverse plane, and to leave at least 12 
hits in the DCH. 

Positive kaon identification is required for the tracks 
forming the candidate <p meson. This is based on dE / dx 
information from the DCH and SVT, and the Cherenkov 
angle and the number of photons measured with the 
DIRC. The kaon selection is based on the likelihood cal- 
culated for each detector component and uses, for each 
track, the ratio of likelihoods for the pion and the kaon 
mass hypotheses, L^/ Lk- If this ratio is less than unity 
for at least one of the detector subsystems, the particle is 
selected as a "loose" kaon candidate. A "tight" identifi- 
cation criterion is also used in the analysis, based on the 
product of the likelihoods for each detector component. 
In this case the track is considered a kaon if the ratio 
of these product likelihoods for the pion- and kaon-mass 
hypotheses is less than unity. 

Three charged tracks originating from a common ver- 
tex are combined to form a Df candidate. Two oppo- 
sitely charged tracks must be identified as kaons with the 
"loose" criterion, and at least one of them must pass the 
"tight" criterion. No identification criteria are applied 
to the pion from Df decay. The reconstructed invariant 
mass of the K + K~ candidates must be within 8MeV/c 2 
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of the nominal <j> mass |J. In the decay Df — > </>7r + , 
the </> meson is polarized longitudinally and therefore the 
angular distribution of the kaons has a cos 2 Or depen- 
dence, where Or is the angle between the K + and Df 
in the <j> rest frame. We require \cos9h\ > 0.3, which 
Monte Carlo studies show retains 97% of the signal while 
rejecting about 30% of the background. 

With these requirements, signals for 4>n + and 

the Cabibbo-suppressed decay D + — > (f)ir + are readily ob- 
served (Fig. ||a) . The Df and D + peaks are both fit with 
single Gaussian distributions with a common free width. 
We model the combinatorial background with an expo- 
nential function. From the fit a D+ signal of 47794 ±311 
events is found with a mass difference m(Df) — m(D + ) 
of 98.4 ± 0.1 ± 0.3 MeV/c 2 . The first error on the latter is 
statistical, and the second is systematic, obtained from a 
study of the mass difference as a function of momentum 
in both data and Monte Carlo simulation. Although the 
uncertainties in the absolute mass scale are on the order 
of several MeV/c 2 , the systematic error in the determina- 
tion of the Z?+ and D + mass difference is much smaller, 
since many sources of error cancel. 

Candidate D* s + mesons are reconstructed in the decay 
D* + — > -D+7, with the subsequent decay Df — + (f>ir + . 
Df candidates are selected by requiring the 4>tt invariant 
mass to be within 2.5 standard deviations (a) of the fit- 
ted peak value. These Df candidates are then combined 
with photon candidates in the event. Photon candidates 
are required to satisfy E 1 > 50 MeV, where E 1 is the pho- 
ton energy in the laboratory frame, and E* > 110 MeV, 
where E* is the photon energy in the T(4S) center-of- 
mass. When combined with any other photon in the 
event, the photon candidate should not form a 7T°, defined 
by a total center-of-mass energy E*^ > 200 MeV and an 
invariant mass 115 < M 77 < 155MeV/c 2 . The distribu- 
tion of the mass difference AM = M(Dfj) - M(D+) is 
shown in Fig. ||b. 

The AM distribution of the signal is parameterized 
with an asymmetric function to account for energy leak- 
age and calorimeter shower shape fluctuations. The sig- 
nal is modeled with a Crystal Ball function j?J , which in- 
corporates a Gaussian core with a power-law tail toward 
lower masses. For the background, a threshold function 

/(AM) =pt(AM -p 2 )P3 e P4(AM- P2 ) 

is used, where the four parameters pi are free in the 
fit. After ensuring that the connection point between 
the Gaussian and power-law tail does not depend on mo- 
mentum and agrees with Monte Carlo simulation, this 
parameter has been fixed to 0.89er in the final fit. A 
signal with 14392 ± 376 D* + events is observed. 



IV. EXTRACTION OF £>r )+ MOMENTUM 
SPECTRA 

The momentum spectrum of mesons in the e + e~ 
center-of-mass frame is extracted by fitting the (jm in- 




1.85 1.9 1.95 2 

M((jm) (GeV/c 2 ) 

N 4000 



> 
33 

"1 3000 - 



1000 - 








0.05 0.1 0.15 0.2 0.25 
M(Dj)-M(D s ) (GeV/c 2 ) 



FIG. 2: (a) The <f>Ti invariant mass spectrum. In addition to 
the Df signal, candidates for the Cabibbo-suppressed decay 
D+ - <h-K + axe also observed. The fit function is a single 
Gaussian for each peak, with widths constrained to be equal, 
plus an exponential background, (b) Distribution of the mass 
difference AM = M D + - M D +. The fit function is a Crys- 
tal Ball function for the signal plus a threshold function, as 
described in the text. 



variant mass distribution for 24 ranges of Df candidate 
momentum. These ranges are 200 MeV/c wide, which is 
much larger than the momentum resolution (w 6 MeV/c) . 
The same function with two single-Gaussians described 
above for the fit to the full mass distribution is used as 
well for the individual momentum bins. Since there are 
many more events in the on-resonance data sample, the 
number of Df in the off-resonance data is extracted with 
the Gaussian parameters (M D + , M D + and er) fixed to the 
values obtained from the on-resonance data. 

The center-of-mass momentum spectrum for D* + 
mesons is extracted by fitting the AM invariant mass 
distribution in 250 MeV/c wide D* + momentum ranges. 
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V. INCLUSIVE BRANCHING FRACTIONS 

The and D* + production cross sections in qq con- 
tinuum are obtained by integrating the momentum spec- 
tra obtained from the off-resonance data. This gives 



a(e + e 



DfX)xB{Dt 



d; ± x)xB(d' s 



= 7.55±0.20±0.34 pb, 



5.8±0.7±0.5 pb, 



where the first error is statistical and second systematic. 
Sources of systematic error are listed in Table |. These 
include the statistical precision of the Monte Carlo de- 
termination of the efficiency, the luminosity uncertainty, 
and contributions from residual uncertainties on track- 
ing (1.2% per track) and particle identification efficien- 
cies, which are determined from control samples in data. 
In addition, for the D* X measurement, there are con- 
tributions from the uncertain signal shape, and residual 
uncertainties on the photon and ir° veto efficiencies, again 
determined with control samples. 

(*)+ 

In order to determine the momentum spectra for Dg 
mesons from B meson decays, the off-resonance data are 
scaled by the on- to off-resonance luminosity ratio and 
then subtracted bin-by-bin from the on-resonance data. 
Integrating the resulting spectrum after continuum sub- 
traction and efficiency correction gives a total yield 
from B meson decays of 87711 ± 1485 events. This cor- 
responds to an inclusive branching fraction of 



B(B -> D+X) 



(10.93±0.19±0.58)x 



(3.6 ±0.9)% 



B{Dl 



%. 



FIG. 3: Efficiency-corrected center-of-mass momentum spec- 
tra for (a) D+ and (b) D* s + for on-resonance (filled circles) 
and scaled off-resonance data (open circles). 



Likewise, the total D* + yield from B meson decays is 
60047 ± 6201 events, leading to the inclusive branching 
fraction of 



B(B - D*+X) = 



(7.9±0.8±0.7)x 



(3.6 ±0.9)% 
B{Df 07T+) 



%. 



We use a larger range because the D* + yield is lower. 
The AM distributions are modeled with a Crystal Ball 
function for the signal and a threshold function for the 
background as described above for the fit to the full dis- 
tribution. The off-resonance data are again fit with the 
Gaussian parameters (x and a) fixed to the values ob- 
tained from the on-resonance data. 



The efficiency e, obtained from Monte Carlo simulation 
of BB and cc events, varies as a function of the D s 
center-of-mass momentum p* . The efficiency ranges from 
20% (5%) when the D+ (D*+) is at rest to 40% (20%) 
for p* — 5 GeV/c. The efficiency-corrected momentum 
spectra of Df and D* + are shown in Fig. 0. 



In the results above, the first error is statistical and 
the second is systematic. The dominant error, due to 
the uncertainty in the Dj — > 4>tt + branching fraction of 
(3.6 ± 0.9)% 0, is shown separately. It is important to 
note that, with this method, the result is independent 
of any assumption regarding the shape of the fragmen- 
tation function. The various contributions to the sys- 
tematic error are listed in Table Q. In addition to the 
sources already noted above, the uncertainty in the shape 
of the background impacts this measurement, particu- 
larly in the lower momentum bins. This contribution is 
estimated with the use of different parameterizations for 
the background shape and different methods for handling 
the continuum subtraction. The efficiency variation over 
the width of the momentum bins is also included as an 
additional systematic error. 
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TABLE I: Systematic errors for cross section and branching 
fraction measurements. 



1YVO 


Fractional 


Error 


(%) 




Continuum 


B Decays 






D*+X 


D+X 


D*+X 


— - 

Signal snape 




3.0 


0.5 


3.0 


Rq f\r cfvcw ~\Yi(\ QI l r>+~ t" £1 ft" l fXTl 






0.4 


4.2 


lY'l i"iyi F o C ar\r\ ctdticfifc 
IVlLUlLfc! V_yd.1 1L) a LdLlo LIL-O 


1.0 


4.8 


2.5 


4.2 


R l n \A7i fi "t" n 

J_J111 WILllill 






1.4 


2.0 


lotal lor JJ S yield 


1.0 


5.7 


2.9 


7.0 


Luminosity /~N(BB) 


1.5 


1.5 


1.6 


1.6 




1.6 


1.6 


1.6 


1.6 


Particle identification 


1.0 


1.0 


1.0 


1.0 


Tracking efficiency 


3.6 


3.6 


3.6 


3.6 


B{D*+ D+ 7 ) 




2.7 




2.7 


Photon efficiency 




1.3 




1.3 


7T° veto 




2.7 




2.7 


Total systematic error 


4.5 


8.2 


5.3 


9.0 



VI. 



FITS TO D 



(*)+ 



MOMENTUM SPECTRA 



(*)+ 

By fitting the D s momentum spectrum, relative 
branching fractions of B decays to different final states 

containing Da mesons are obtained. In the T(AS) rest 

(*)+ 

frame, two-body B decays produce Ds mesons with 
a momentum spectrum about 300 MeV/c wide. In B de- 

(*)+ 

cays, the D s momentum spectrum is essentially gov- 

(*)+ — 

erned by the production of direct D s . Other cs states 
(with L = 1), such as L»+ 1 (2536) and D*^(2573), pri- 
marily decay to D^K. Because D*+ decays to D +r j or 
D+~ir°, the D+ momentum distribution is slightly broader 
and shifted downward compared to direct production 
from B -> D+X. 

Three different sources of D s *^ + mesons in B decays 
are considered for the fits to the momentum spectra. 

(1) B — * Dy' + D^ decays. The relative branching 
fractions of the individual channels can be taken 
either from existing measurements H or from pre- 
dictions that assume factorization || [h], [ll[] . The 
fit is performed for both cases, with the assump- 
tion f D *+ = f D + for the theoretical models, where 

f D (')+ are the D s decay constants. 

(2) B -> D S *^D** decays_ Four Upstates are con- 
sidered: Dl{j = 1/2), L>i(2420), = 1/2) and 
^(2460). Observation of B -> D { s * )+ D** decays 
was recently reported by CLEO |Q. 

(3) Three-body B -> D ( s * )+ D^n/p/to decays. Since 
little is known about these decays, they are at- 
tributed equal weights, and the momentum distri- 
butions are generated according to phase space. 

Minimum-x fits to the D s momentum spectra are 
performed, where the total number of D s events and 
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FIG. 4: Fit results for (a) D+ and (b) D^+ momentum spec- 
tra. The data are dots with error bars, and the histograms 
are the components of the fit function described in the text. 
Type (1) is B -> Di* )+ D M , Type (2) is B -> D^+D" and 
Type (3) is B -> D ( s * )+ D^tt/ p/uj. The solid histogram is the 
sum of the three components. 



the fractions of the source (1) and (2) contributions are 
free parameters. From the fits to the D+ and D*+ spec- 
tra, the ratios of two-body modes [source (1)] to the total 
inclusive rate arc determined to be 



EB{B -> D 



B(B 
SB{B- 



-» D+X) 



B(B -> D* s +X) 



= (46.4 ±1.3 ±1.4 ±0.6)%, 
= (53.3 ±3.7 ±3.1 ±2.1)%. 



The first error is statistical. The second error repre- 
sents the systematic error due to the limited Monte Carlo 
statistics and the background parameterization. 

The last error is due to the model uncertainty. It is 
obtained by varying the relative fractions of the modes 
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(*) + 

contributing to each source of Dg listed above. The 
fit is performed with alternative assumptions for the 
relative contributions of the modes in source (1) taken 
from theoretical predictions and measurements. Dif- 
ferent weights for B -> D+D** and B^D* S +D**, as 
well as different relative branching fractions of the four 
modes within source (2), are used. For source (3), ei- 
ther B -> D^D^ir or B -> D { s * } 5W p / 'uj is assumed 
to be dominant. The x 2 °f the fit for the inclusive 
D* + momentum spectrum is lowest when the contri- 
bution of B -> D ( * ] D^p/ lu is dominant compared to 

B -► D { * ] D^\. The results of the fits to the D { s * )+ 
momentum spectra are shown in Fig. |] under the as- 
sumption of equal weights for the individual contribu- 
tions within sources (2) and (3), and with the weights of 
the individual modes of source (1) taken from 0|. 

The sum of branching fractions for the two-body B — > 
D ( S * ] D^ decays are obtained from the fits to the D\. ' 
momentum spectra, where the yield from each source is 
a free parameter. We find 

EB(B -> L>(* )+ Z3 ( * ) ) = (5.07 ± 0.14 ± 0.30 ± 1.27)%, 
SB(B -> D*+ Z? w ) = (4.1 ±0.2 ±0.4 ±1.0)%, 

where the first error is statistical, the second is system- 
atic, and the third is due to the Df — * 4>n + branching 
fraction uncertainty. The systematic error includes con- 

(*) + 

tributions from the B — > D s X branching fractions, 
the relative contributions of source (1), and model de- 
pendence of the source spectra. The sum of the two-body 
modes is reasonably separated in the momentum spectra 
from the other components. Therefore, the fractional er- 
ror on the sum of the two-body modes is smaller than the 



fractional error on the B — > D s X branching fraction 
or the relative two-body branching ratio. 



VII. SUMMARY 

In summary, the branching fractions for inclusive B — > 

(*)+ 

D s ' X production have been determined as well as the 

(*)+ 

D s production cross sections from continuum events 
at center-of-mass energies about 40 MeV below the T(4S) 
mass. Our more precise results for the Df are in agree- 
ment with previous measurements ^ [b| , while the D* + 
measurements are new. In contrast to previous results, 
our measurements do not rely on any assumptions re- 
garding the shape of the fragmentation function. Fi- 
(*)+ 

nally, fits to the D s momentum spectra provide rel- 
ative yields and branching fractions for two-body B — » 
£)(*)+£(*) an d B _> £>*+/J(*) decays. The mass differ- 
ence m(D+) — m(D + ) has also been measured. 
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